A period distribution of X-ray binaries observed in the central 
region of M31 with Chandra and HST 



R. Barnard, J. L. Galache, and M. R. Garcia 

Harvard- Smithsonian Center for Astrophysics, 60 Garden St, Cambridge MA 02138 
O ■ and 



(N 

o 

(N 



m 



s 



N. Nooraee, and P. J. Callanan 
Dublin Institute for Advanced Studies, Dublin, Ireland; University of Cork, Cork, Ireland 



W . A. Zezas 

HH \ University of Crete, Greece; Harvard-Smithsonian Center for Astrophysics 

r^ ' and 

^1 S. S. Murray 

f/3 ■ Johns Hopkins University, Baltimore, Maryland; Harvard- Smithsonian Center for Astrophysics 

m 

Q^ : ABSTRACT 

^^ . Almost all Galactic black hole binaries with low mass donor stars are transient X-ray sources; 

we expect most of the X-ray transients observed in external galaxies to be black hole binaries also. 
Obtaining period estimates for extra-galactic transients is challenging, but the resulting period 
distribution is an important tool for modeling the evolution history of the host galaxy. We have 



m 

fvj . obtained periods, or upper limits, for 12 transients in M31, using an updated relation between 

the optical and X-ray luminosities. We have monitored the central region of M31 with Chandra 
for the last '^12 years, and followed up promising transients with HST; Aa B magnitude limits 
for optical counterparts are 26-29, depending on crowding. We obtain period estimates for each 

^^ , transient for both neutron star and black hole accretors. Periods range from <0.4 to 490±90 

hours (<0.97 to <175 hrs if all are BH systems). These M31 transients appear to be somewhat 
skewed towards shorter periods than the Milky Way (MW) transients; indeed, comparing the 
M31 and MW transients with survival analysis techniques used to account for some data with 
only upper limits yield probabilities of ~0. 02-0. 08 that the two populations are drawn from the 
same distribution. We also checked for a correlation between orbital period and distance from 
the nucleus, finding a 12% probability of no correlation. Further observations of M31 transients 
will strengthen these results. 

Subject headings: x-rays: general — x-rays: binaries — black hole physics 

1. Introduction of that galaxy. One way of characterizing such 

populations is the distribution of their orbital pe- 
X-ray binaries (XBs) represent the end-points ^jo^s; however, the estimation orbital periods for 

of stellar evolution, and the XB populatfon of a extra-galactic XBs is extremely challenging, unless 

galaxy is informed by the star formation history 



they happen to be at high enough incHnations to 
exhibit dips and/or echpses. 

The most accessible extra-galactic population 
for period studies is the M31 X-ray transient 
group. These are expected to be low mass X-ray 
binaries (LMXBs) with black hole (BH) primaries, 
since t he Galactic BH LMXBs are nea rly all tran- 
sient ( McClintock fc RemillardI 120061 ): transient 
neutron star LMXBs also exist, but the major- 
ity of transients contain BH primaries (see e.g. 
Chen et al.lll997() . X-ray and optical observations 
of these transients allow us to estimate their size, 
and therefore period, using the methods described 
below. Comparing the period distributions of the 
M31 and Galactic BH transient LMXBs may shed 
light on possible differences in the star formation 
and binary evolution histories of these galaxies. 

Ivan Paradijs fc McClintockl (|l994l hereafter 
VPM94) found the absolute V band magnitudes 
of 18 Galactic LMXBs to vary between —5 and 
+5, depending on the size of the accretion disk, 
and X-ray luminosity. To do this, they required a 
sample with well known distances, and therefore 
limited the sample to systems exhibiting X-ray 
bursts with radius expansion, soft X-ray transients 
and Z-sources. The Eddington limit was assumed 
to be 2.5x10^* erg s~^. They found that the op- 
tical counterparts to bursters were fainter than 
those of non-bursters; this was a first clue that 
the optical and X-ray luminosities were linked, 
because bursters tend to have low X-ray lumi- 
nosities, while the non-bursting Z-sources operate 
near the Eddington limit. Furthermore, previ- 
ous work suggested that My was proportional 



to the orbital period (e.g. ICheva 



Nieto et al.lll990l: lAnderson et al.lll993l ): this led 



ier et al.l [l985[ 



VPM94 to postulate that the relation between 
X-ray and optical luminosities (Lx and Ly re- 
spectively) was caused by X-ray reprocessing in 
the accretion disk. 

The VPM94 toy model suggested that Ly oc 

1/2 

L-^ a, where a is the binary separation. One 
might expect Ly to be proportional to a^, since 
that is related to the area of the disk assuming 
that the disc radius is proportional to a. How- 
ever, the disk temperature decreases with radius: 
T" ex Lx/a^, and Ly oc a^ T" where a ~2 for the 
observed temperature range in LMXBs (VPM94) . 

Kepler's third law gives a oc P^/^ (Mi -I- A'h)^^^; 



VPM94 simplified this by assuming that the 
masses of the two stars were similar, and defined 
S as (Lx/iEdd) (-P/1 hour) ' . Their empirical 
relation 

My = 1.57 (±0.24) - 2.27 (±0.32) log S (1) 

provided a good match to observations over 3 or- 
ders of magnitude in S. However, their sample in- 
cludes black hole and neutron star accretors, and 
a mix of different X-ray spectral states. 

More recently iRussell et al.l (|2006l ) examined 
the correlations between optical/infra-red and X- 
ray luminosities for neutron star and black hole 
X-ray binaries separately, also factoring in the 
masses of the accretor and donor stars (i.e. Mi 
and M2). They found that the broadband BVRI 
luminosity (Loir) correlated well with L^^ over 
8 orders of magnitude in Lx for black hole bina- 
ries, and also neutron star binaries in the hard 
state; the black hole normalization was ^10 times 
higher than the neutron star normalization, con- 
sistent with the expectation that black hole sys- 
tems have larger disks for a given orbital period 
because {Mi + M2) is larger for black holes than 
for neutron stars. 

Moss et al. (2012, in prep, hereafter M12) fo- 
cused on the My vs. E relation for black hole 
binaries in the high state at the peak of their out- 
bursts, with the aim of providing the best period 
estimates for X-ray transients. They used X-ray 
and optical data for transients near their outburst 
maximum. Their new empirical relation 

My = 0.84 (±0.30) - 2.36 (±0.30) log E (2) 

tends to result in a shorter period than the relation 
described in Equation 1. We note that the only 
significant difference between Equations 1 and 2 
is in the normalization, again caused by the larger 
disks in the black hole systems. 

We have conducted a 12 year campaign to mon- 
itor the central region of M31 for bright X-ray 
transients. Over 130 Chandra observations were 
made over this time, averaging one per month; 
no observations are made when M31 is too close 
to the sun. Bright new transients are followed 
up with HST ACS WFC observations using the 
F435W filter that approximates the B band. One 
observation was made shortly after the transient 
outburst, and another observation was made ~6 



months later, if there were no sere ndipitous cover 
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age in an carhcr obse rvation (see 
[2005a, c, 2006; Barn ard et al 
et al. 2012). In this work, we present our analy- 
sis of 12 X-ray transients observed with Chandra 
ACIS and HST. Herein we present the discovery 
of 2 new optical counterparts, and provide revised 
period estimates for all 12 systems. 

X-ray transient behaviour appears to be gov- 
erned by thermal-viscous instability in the ac- 
cretion disc, with hot and cold stable phases, 
along with an unstable transitional phase (see e.g. 
Lasotal l200ll and references within). The out- 



bursts occur when the accretion disc enters the hot 
stable phase; however, the resulting X-ray illumi- 
nation prevents the disc from returning to the cold 
state until the X-ray l uminosity is sufficiently di- 
minished (|King fc Ritter 1998. hereafter KR98). 
KR98 showed that we may expect exponential de- 
cay lightcurves so long as the the whole disc is 
irradiated with X-rays; otherwise, linear decay is 
expected. For long transient outbursts, KR98 ex- 
pect the decay to switch from expo nential to linear 
after a few e- folding times. Shahbaz et al.l ( 19981 ) 
analyzed the decay lightcurves of various tran- 
sients, and found good agreement with the theory 
of KR98; furthermore, they found that systems 
with multiple outbursts tended to exhibit the same 
type of decay during each outburst. 

2. Observations and data reduction 

We provide a journal of HST observations in 
Table [H For each HST ACS observation we pro- 
vide a number, the observation ID, the date, the 
exposure time, the pointing, and the target (la- 
beled T1-T12 in order of increasing RA); we also 
provide the closest ACIS observation, along with 
the time of observation with respect to the HST 
observation. If the transient is in a region that 
we have already observed, then a dedicated "off" 
observation is unnecessary. 

2.1. X-ray analysis 

We have examined 90 Chandra ACIS observa- 
tions of the central region of M31, spaced over 
~12 years; these observations comprise a monitor- 
ing programme that is designed to find transient 
X-ray sources in the central region of M31. The 
roll-angle of each observation was unconstrained. 



resulting in an approximately circular field of view 
with 20' radius. We determined the position of 
each source from a merged 0.3-7 keV ACIS im- 
age, using the iraf tool imcentroid. 

The merged ACIS image was registered to the 
Field 5 B band image from the Local Group 



Gala xy Survey observations of M31 (jMassev et al. 



20061 ). using 27 X-ray bright globular clusters. We 
were unable to register the Chandra positions with 
the HST observations directly because only a few 
of the HST fields contained X-ray bright globular 
clusters. After refining the positions of the globu- 
lar cluster in the Chandra and LGS B band images 
with IMCENTROID, we found the sky co-ordinates 
of the new positions with xy2sky, then registered 
the X-ray positions to the LGS B band image us- 
ing the IRAF task CCMAP. The mean r.m.s. un- 
certainty in position due to registration was 0.11" 
in R.A. and 0.09" in Dec. 

We extracted source and background lightcurves 
and spectra, along with the appropriate response 
files, for each observation of all 12 sources. Source 
and background regions were circular, with the 
same radius; the radius was optimized for crowd- 
ing and off-axis angle, and ranged over ^2-30". 
The location of each background region was cho- 
sen to be near to the source region, but not so close 
that it might be contaminated by source flux. 

We created long term, calibrated lightcurves for 
each source, obtaining a 0.3-10 keV unabsorbed 
luminosity for each observation. We used these 
lightcurves to determine the overall variability of 
each X-ray source, and the X-ray luminosity at the 
times closest to each HST observation. Since we 
expect our transients to be black hole X-ray bina- 
ries in the high/soft state, we initially estimated 
the conversion from flux to luminosity by assum- 
ing a disk blackbody emission spectrum with in- 
ner disk temperature kTin = 1 keV, and A^h = 
6.7x10^° atom cm^^ (the Galactic line of sight 
column density), then determining the unabsorbed 
0.3-10 keV luminosity equivalent to 1 count s~^. 
After correcting for the exposure, vignetting and 
background, multiplying the source intensity by 
this conversion factor gave the source luminosity. 
Source spectra with >200 net counts were freely 
fltted with disk blackbody models. When good 
spectral fits were found for a particular observa- 
tion of a source, the parameters of these fits re- 
placed the default parameters when estimating the 



source luminosity for observations of the source 
with <200 photons. 

2.2. Optical analysis 

For each ACS observation, we examined two 
types of images provided by the STSCI. Each ob- 
servation had contained one drizzled (DRZ) image 
and 4-8 flat-fielded (FLT) images. The drizzled 
image is the final product: the flat fielded images 
are combined, the sky background is subtracted, 
and cosmic rays are removed; the image is scaled 
by intensity (electron s~^). The FLT images were 
corrected for instrumental effects (dark, bias, and 
flat); however, the sky was not subtracted, and 
cosmic rays were not removed; furthermore, the 
image was scaled to the total number of electrons 
received over the observation. We used the DRZ 
images to create the difference images described 
below, but the FLT files for photometry since 
DAOPHOT requires an accurate measurement of 
the sky background. 

Each drizzled HST observation was registered 
to the LGS B band image, using the method de- 
scribed in the previous section. We then made a 
difference image for each pair of on/off observa- 
tions. To achieve this, we first used the iraf task 
WREGISTER to ensure that the pixels in each image 
had the same co-ordinate system and orientation. 
We then used the farith ftool to create the dif- 
ference image. We examined the difference image 
for each observed variable X-ray source, looking 
for significant changes within a 3cr radius, includ- 
ing the uncertainty in the centroid position, and 
the r.m.s uncertainty in registration of the Chan- 
dra image to the LGS B band image; the r.m.s 
uncertainties in registering the HST images to the 
LGS image were negligible. 

If no counterpart was found in the difference 
image, then we estimated the Aa upper limit to 
the B magnitude of an undetected star. To do 
this, we summed the counts in each of the FLT 
images, within a circle with 2 pixel radius at the 
position of the X-ray source, for a total of Ctot 
counts over T seconds. The maximum intensity 
of the counterpart is then 4(Ctot) /T, which we 
convert to Vega B magnitude via 

B = -2.5 log [4 (Ctot)"-^ /t] + 25.76. (3) 

In order to ensure that our B magnitude limits 



were believable, we ensured that DAOPHOT was 
able to detect several stars with similar or fainter 
magnitudes in the FLT images. For images where 
the sky level varied significantly, we restricted the 
search for similarly faint stars to areas with similar 
sky levels as the target. 

Unfortunately, we did not have the simulta- 
neous B VRI observat i ons to utilize the relation 
found bv lRussell et al.l pOOGl ): however, we did ob- 
tain period estimates using the VPM94 and M12 
methods. Both the VPM94 and M12 methods 
for estimating the orbital period require the abso- 
lute magnitude in the V band. We assumed that 
our F435W flux represented the B band flux, and 
used the relation B— V = 0.09 that was empiri- 
cally derived for LMXBs (JLiu et al.l [20011) . Sys- 



tematic studies of X-ray scattering halos around 
Galactic LMXBs have shown that Ay ~ JVh / 
1.8 X 10^1 atom cm^^ (jPredehl fc Schmitd 119951) 



while Av/E{B - V) ~3 (see e.g. iCardelli et al. 
1989); hence Ab ~ ^h x (l-hl/3)/1.8 x lO^i. 
Fi nally, the distance mod ulus for M31 is 24.47 
(jStanek fc GarnavichllT998h . Hence 



My = B-t-0.09-iVHx(l -I- 1/3) /1.8x 10^1-24.47. 

(4) 

2.3. Estimating the periods 

When calculating the periods derived from the 
M12 relation, we were able to include uncertain- 
ties in the relation. We achieved this by randomly 
altering the My of every Galactic black hole bi- 
nary high state in the Moss et al (2012) sample ten 
thousand times, drawn from a normal distribution 
scaled by the uncertainties, and obtaining a new 
My vs. log E relation for each iteration; these My 
distributions had means and spreads that were 
consistent with the data. We then obtained the 
log S corresponding to the My of each target in 
every iteration, and derived the mean log E and 
uncertainties for each target by fitting a Gaussian 
to the resulting log E distribution. Unfortunately, 
we did not have access to sufficient information to 
include uncertainties in the VPM94 relation; in- 
stead, we propagated the uncertainties in My and 
Lx assuming that the relation was fixed. 

3. Results 

Figure [1] shows the locations of the 12 transients 
on a log-scaled B-band image of the central re- 



gion of M31, provided by the LGS (Massey et al., 
2006). Table [2] gives the positions of each target, 
and summarizes the results discussed later in this 
section; we provide the 0.3-10 keV luminosity dur- 
ing the closest ACIS observation to the HST on 
observation, our estimated 0.3-10 keV luminosity 
at the time of the HST on observation, the B mag- 
nitude (or 4 a upper limit) for the counterpart, 
and period estimates using the empirical relations 
of VPM94 and M12. We present the ~12 year, 
0.3-10 keV unabsorbed luminosity lightcurves of 
the 12 transients in Fig. [51 

3.1. Analysis of the Chandra observations 

Only four of the targets had >200 photons in 
the Chandra observation that was closest to the 
HST on exposure: Tl, T2, T5, and T9. Spec- 
tral fit parameters are provided in Table [3] for ab- 
sorbed disk blackbody models. The T9 spectrum 
was piled up; we account for this in Noorae et 
al. (2012, submitted). Only Tl apparently ex- 
perienced greater than Galactic absorption, with 
iVn = 2.5±1.9xl0^^ atom cm~^ (uncertainties are 
quoted at the 90% confidence level). We find that 
Ab = 1.9±1.4 for Tl, and 0.50 for the remaining 
transients. 

VPM94 obtained X-ray luminosities from the 
Ariel V X-ray c atalog, and Ariel V e mployed the 
2-10 keV band ([Warwick et al.lll981l ). Hence we 
estimated the 2-10 keV luminosity of each target 
at the time of the HST observation. For those 
transients with freely fitted spectra, we used the 
i2-io/^o.3-io factors provided in Table [31 For 
transients with no freely fitted spectra, we assume 
kTin = 1 keV, giving ^2-10/^0.3-10 = 0.48. 

However, it is possible that certain transients 
were in their low/hard state at the time of the 
HST on observation; the transition from low state 
to high state during the rise of the outb urst occurs 



at a luminosity 5,10% E ddington (Gladstone et al 



2OO7I : [Tang et al.ll2011[ ). while the transition from 
high state to low state during the decay occurs 
at ^^2% Edding ton due to hysteresis (see e.g. 
Maccarondl2003l . and references within) . The low 



state may be characterized by a power law with 
photon index F '^ 1.7; assuming a power law with 
F = 1.7 rather than a disk blackbody with kTjn 
= 1 keV results in 2-10 keV luminosities that are 
~75% higher, and periods that are ^^35% shorter. 
Since 8 of the 12 transients only yielded upper 



limits to their periods, we present the more con- 
servative upper limits obtained by assuming high 
state spectra. 

3.2. Analysis of the HST observations 

We present the difference images for T1-T12 
in Fig. [31 showing the difference in intensity be- 
tween the drizzled images of the on and off obser- 
vations; white is positive, and black is negative. 
The 1(7 and 3(T uncertainty ellipses are shown, in- 
cluding uncertainties in X-ray position, and X-ray 
registration to the LGS B image. North is up, 
east is left. Each image shows a 5"x5" field. The 
background level varied considerably, with the nu- 
clear region being drastically more crowded than 
the outer regions; as a result, 4it detection limits 
ranged over B^26-29. 

3.3. Discussion of individual sources 

3.3.1. Tl 

Chandra, Swift, HST and Einste in observations 
of Tl are discussed in detail by [Barnard et al 



([2OIID . Tl has no ACS off image with the F435W 
filter, because the WFC was not in use at the 
time. Instead, we used a serendipitous 81 s V 
band ACS/WFC observation (ID J92ga5naq, 2005- 
02-17, PI Crotts). There are many black stars in 
the difference image, because they are brighter in 
the V band than in the B band; however, there 
are no white stars to signify the transient, and we 
estimate that B >28.7. The sky background was 
extremely low, and we detected several stars over 
the image with B >28.7 at the 4cr level in one FLT 
image. Since there was no B band off image, we 
examined the 8 individual FLT images to deter- 
mine whether there were any detections of stars 
within the 3ct ellipse; this allowed us to discrimi- 
nate between stars and e.g. cosmic rays. We found 
no stellar candidates that were detected in >3 out 
of 8 FLT images, and conclude that there were no 
bright B band stars within the 3cr ellipse. 

For Tl, the closest ACIS observation to the 
HST on image was Obs 7139 (day -2500 in Fig. 
2), taken 27 days previously; however, Tl was 
not covered by any ACIS observations for >100 
days either side of the 2006 outburst. In a later 
outburst (day ~3900), the 0.3-10 keV luminosity 
went from 8.6x10^'^ to 0.4x10^^ erg s"^ over 27 
days. Assuming exponential decay, the e-folding 



time was 8.8 days. We applied the same exponen- 
tial decay to the 2006 outburst, and found a 0.3-10 
keV luminosity of 0.5x10'^^ erg s~^ at the time of 
the HST on observation. For the peak spectrum, 
we obtained a best fit kTi„ of 1.5±0.3 keV; as a 
result, L2-io/-^o.3-io = 0.55. 

3.3.2. T2 

The T2 outburst (day ~2550) lasted ~150 days. 
The X-ray peak was during Obs 7140, 32 days 
prior to the HST on observation. Fitting the 
lightcurve with exponential decay yielded an e- 
folding time of 38 days; however, the residu- 
als were large, and the x^/dof was unacceptable 
(47/3). We also fitted the lightcurve with power 
law and straight line models; the best fit values 
had worse x^/dof than the exponential fit. There 
was insufficient data to estimate the intensity for 
a model where the decay went from exponential 
to linear, hence the exponential decay model pro- 
vided the best estimate of the 0.3-10 keV luminos- 
ity of T2 at the time of the HST on observation. 
Using the best fit e-folding time of ^^40 days, we 
obtained a 0.3-10 keV luminosity of ~6xl0^^ erg 
s~^ for the HST on observation. The best fit disk 
blackbody model for the peak spectrum yielded 
kTin = 0.54 ±0.04 keV; this results in a 2-10 kcV 
luminosity of just 1.5x10^^^ erg s~^. 

The positioning of the T2 on observation was 
off by >1', due to the misidentification of the guide 
star. As a result, it was impossible to accurately 
register the FLT images to the drizzled image, and 
this limited our photometric accuracy. The differ- 
ence image shows one particularly variable star; 
its B magnitude changed from ^24.7 to ^^25. 6 be- 
tween the on and off observations. We therefore 
consider this star to be the most likely counter- 
part. Since the our optical counterpart to T2 had a 
B magnitude ~25.6 even after the X-rays had dis- 
appeared, we removed this component when cal- 
culating the period. 

3.3.3. T3 

T3 was first studied bv lWilliams et al.l (|2005bl ). 



who found no counterpart and concluded that B 
>25.5. However, our method of estimating the 
magnitude limit for the counterpart yielded B 
>27.2. 

T3 exhibited two outbursts, ^2000 days apart; 



the HST observations were triggered by the first 
outburst. This outburst had an observed peak 
0.3-10 keV luminosity of 3.7±0.3xl037 erg s-\ 22 
days before the HST on observation (day ~1700); 
however, the next observation of T3 was ~150 
days afterward, where the 0.3-10 keV luminos- 
ity had fallen to 0.35±0.09 xlO^^ erg 8"^ If 
we assume that these observations cover a single 
exponential decay, then the e-folding time is 68 
days. The 0.3-10 keV luminosity would then be 
-2.73x103"^ erg s'^ at the time of the HST on 
observation. 

3.3.4. T4 

T4 was mostly quiescent, but exhibited a ~1500 
day period of X-ray activity, starting on day 
~1500 where the intensity varied by a factor ~20. 
The HST on observation was triggered at day 
-2500. The 0.3-10 keV luminosity of T4 was 
>1.5xl0^^ erg s~^ for the Chandra observations 
closest to the HST on observation, and ~0.4xl0'^^ 
erg s~^ at the time of the off observation. 

We found no evidence for stellar variability in 
the difference image of T4; we set B >28.9. 

3.3.5. T5 

T5 was first discussed bv lWilliams et al.l (|2006( ) , 
who found no counterpart, and set a magnitude 
limit of B > 24.7 due to the high sky background 
caused by proximity to the nuclear region. Our 
careful analysis of the FLT images from this obser- 
vation with DAOPHOT improved the magnitude 
limit to B > 26.6. 

T5 exhibited a single outburst that was ob- 
served for —80 days, starting at day —1750. Ex- 
ponential fits to the decay lightcurve were rejected 
(x^/dof >63/3); however the lightcurve was well 
described by linear decay with slope —0.06 (x^/dof 
— 3.5/3). The HST on observation occurred 29 
days after the observed peak, hence we estimate 
the 0.3-10 keV luminosity to be 4.0x10^^ erg s~^. 
The brightest observation yielded a spectrum with 
kTin = 1.1±0.2, yielding a 2-10 keV luminosity of 
2.5x10^^ ergs-^ 

3.3.6. T6 

We discovered from our long-term ACIS lightcurve 
of T6 that it is not a transient, but does vary by 
a factor >10. Fortunately, T6 was observed 6 



times in our survey, and we were able to make a 
difference image from the two observations with 
the highest contrast in X-ray himinosity (Obser- 
vations 3 and 4 in Table 1). We found no evidence 
for variability in the difference image within the 
3cr ellipse; instead, we estimate B to be >26.3. 

We estimate that the 0.3-10 keV luminosity of 
T6 varied by a factor ~4 between HST obser- 
vations 3 and 4 (as described in Table 1), from 
-0.9x103'^ to -0.2x10^'^ erg s^^ 

3.3.7. Tl 

T7 exhibited a slow, faint outburst that lasted 
~250 days, starting at day ^-^2300. The shape of 
the outburst was unlike the classical fast rise fol- 
lowed by exponential decay (FRED); instead, the 
lightcurve was approximately flat for the four ob- 
servations either side of the HST on image, with 
a 0.3-10 keV luminosity of ~ 0.6x10^'^ erg s~^ 

Since the 0.3-10 keV luminosity of T7 did not 
significantly vary between its on and off observa- 
tions, we cannot expect to detect optical variabil- 
ity in the counterpart. As a result, we must set 
our magnitude limit to that of the brightest star 
within the 3cr position uncertainty ellipse, giving 
B >25.0. 

3.3.8. T8 

T8 was initially quiescent, turned on around 
day 1400, and has been active ever since. The 
0.3-10 keV luminosity was ^0.4x10^^^ erg s~^ at 
the time of the HST on observation, and also at 
the the time of the off observation. 

As with T7, there was no significant difference 
between the X-ray luminosity of T8 at the times of 
the HST on and off observations. Hence we set our 
magnitude limit to be fainter than the brightest 
star within 3ct of the X-ray position, i.e. B >23.9. 

3.3.9. T9 

T9 was brightest X-ray sour ce ever seen in M31 , 
and is discussed in detail in iKaur et al.l (J2012I ) 
and Noorae et al. (2012). The difference image 
shows a clear counterpart within the la region. B 
= 23.91±0.08. No counterpart was observed in the 
"off" image, and B >27.8 at the 4(t level; hence, 
we assume that all of the optical light comes from 
reprocessed X-rays. 



The closest ACTS observation was 14 days after 
the HST on image. However, the observed peak 
occurred in a HRC observation, 34 days before 
the HST on observation. The decay lightcurve ap- 
pears to transition from exponential decay to lin- 
ear decay (see Nooraee et al., 2012); we estimate 
the 0.3-10 keV luminosity to have been ^50x 10'^^ 
erg s~^ at the time of the HST on observation. 
Using a disk blackbody emission model, and cor- 
recting for pile-up effects, yielded kT '^0.8 keV; as 
a result, we estimate the 2-10 keV luminosity to 
have been '^18.6x10'^^ erg s~^ at the time of the 
HST on observation. 

3.3.10. TIO 

TIP is discussed in detail by IWilliams et al 



(|2005al ). who used DAOPHOT and ALLSTAR to 
disentangle the counterpart from its bright neigh- 
bor. The counterpart dominated the blend for the 
on observation, but was significantly fainter in the 
off observation, allowing lWilliams et al.l ( 2005al ) to 
obtain B = 24.90 for the other star. Assuming 
the magnitude of the other star was 24.90 dur- 
ing the on observation yielded B = 24.52±0.08 for 
the counterpart during the on observation, and 
24.95±0.08 during the off observation. As with 
T2, we subtracted the "off" component of the op- 
tical emission when calculating the orbital period 
for TIO. 

The TIO outburst was captured in 5 ACIS 
observations over '^180 days; the shape is flat- 
topped, rather than the typical fast rise followed 
by exponential decay. The closest ACIS obser- 
vation occurred 1 day before the HST on obser- 
vation, with 0.3-10 keV luminosity 3.0±0.2xl03'^ 
erg s~^; given the slow evolution of the outburst, 
we assume the same luminosity at the time of the 
HST on observation. 

3.3.11. Til 

The Til outburst was significantly detected in 
5 ACIS observations over 54 days. The decay was 
well described by an exponential with an e-folding 
time of 33 days (x^/dof = 2/5). We calculated the 
0.3-10 keV luminosity to be ~0. 34x10^^ erg s"^ 

The Til difference image revealed a likely coun- 
terpart within the la ellipse. We found that B 
=24.87±0.09 for the on observation. We found no 
trace of a counterpart in the "off" observation; B 



>28.2 at the Aa level. Hence, we assume that all 
of the optical light comes from, reprocessed X-rays. 

3.3.12. T12 



T12 was first discussed by IWilliams et al 



(|2005bl ). who identify a counterpart despite the 
fact that it is brighter in the off observation than 
in the on observation; it appears as a black star 
in Fig. |3l We disagree with this interpretation, 
since there is no scenario where the disc brightens 
during decay, and instead impose a 4cr upper limit 
of B > 27.9. 

T12 exhibited an apparently double-peaked 
outburst over ^200 days, with the second peak 
brighter than the first. However, the previous ob- 
servation to include T12 was ~150 days earlier, 
meaning that the true peak was not observed. 
The 0.3-10 keV luminosity fell from 2.4±0.3xl03'^ 
to 0.26±0.09xl0'^^ erg s~^ over 61.5 days; if we 
assume exponential decay, then the 0.3-10 keV 
luminosity was ~0.9xl0'^^ erg s~^ at the time of 
the HST on observation. 

3.4. Estimating the periods for the 12 tar- 
gets 

Since the e-folding time for transient decay is 
~2.2 times longer in the optical than in the X-ray 
lightcurves, the VPM94 and M12 relations only 
truly apply to the peak of the outburst, and get 
more inaccurate as the burst evolves. We made 
two sets of period calculations for each source: one 
set estimating the X-ray luminosity of each tran- 
sient at the time of the HST observation, and one 
set estimating the optical luminosity at the time of 
the X-ray peak. Comparison of the two sets of pe- 
riods indicated the robustness of these estimates. 

For Tl, the NS periods agreed within to within 
30%, and the BH periods varied by 40%; the e- 
folding time was particularly short for this system. 
For T9 there was a ~30% discrepancy between the 
periods. For T12, the periods differed by a fac- 
tor ^3; this large discrepancy may be due to the 
double-peaked nature of the T12 outburst. The 
two sets of periods agreed to within ~10% for all 
targets not previously specified. These results sug- 
gest that the observed peaks were generally close 
enough to the true peaks for the VPM94 and M12 
relations to hold. 

The periods listed in Table[2]are those obtained 



by estimating the X-ray luminosity at the time of 
the HST on observation; we provide lu uncertain- 
ties for systems with identified counterparts, and 
3ct upper limits to the periods for systems with 
magnitude limits. We note that some of the tran- 
sients are rather X-ray faint, suggesting that they 
are more likely to be in the low/hard state, rather 
than the thermal high state; the bolometric lumi- 
nosities of these sources would then be somewhat 
higher, leading to shorter periods. However, most 
of these systems only have upper limits to the pe- 
riod already; Til is the only X-ray faint source 
with a counterpart. 



We note that IWilhams et all (|2005bl ) fitted the 
ACIS spectra of TIO despite their low quality 
(^^120 net source counts). They favor power law 
emission models with Nh ^2x10^^ atom cm~^; 
they dismissed disk blackbody emission models be- 
cause the A'h was below the line-of-sight value in- 
stead of fixing A^H- As a result, they found Lx = 
6x10^^ erg s-\ and My = -1.3±0.5, both a fac- 
tor ~2 brighter than our values (see Table [2]); these 
results give a ~130 hr period, considerably longer 
than the ~50 hr period that we obtained. How- 
ever, we checked these spectra and found Ah to be 
consistent with the Galactic line-of-sight absorp- 
tion, and we prefer the results from our assumed 
spectral model. 

3.4.1. Comparison of the M31 and MW black 
hole transient periods 

In Fig. |4] we order the transients by increas- 
ing BH period; the y axis serves to separate 
them, but does not reflect any physical param- 
eter. Circles and stars represent NS and BH 
periods respectively, for systems with counter- 
parts; otherwise, we give the 3cr upper limits, 
represented by arrows. We show the period dis- 
tri bution of 15 Milky Way BH binaries listed 



McClintock fc RemillardI (|2006[ ) for compari 



son, represented by squares; the three persistently 
bright BH HMXBs are excluded. Our period dis- 
tribution for the M31 transients is similar to that 
of the Galactic BH LMXBs; however the M31 
transients may tend towards shorter periods. 

We have not included the period for T8 since 
its X-ray magnitude did not vary between the 
HST on and off observations, meaning that our B 
band magnitude limit was set by the brightest star 
within the 3(7 positional uncertainty. This limit is 



brighter than most of the observed counterparts, 
and the resuhing period hmit is not usefuL While 
the X-ray luminosity of T7 did not vary between 
the HST on and off observations either, the B mag- 
nitude limit was fainter than any of the observed 
counterparts so we have included T7. 

We tested the possibility that the periods of 
the transients in M31 are consistent with those in 
the Milky Way using the tests for univariate data 
implemented in the ASURV package; this method 
allows us to compare a censored population, con- 
taining some systems with upper limits only, with 
an uncensored population, where all t he values are 
known, within unc ertainties (see e.g. IPeto fc Petd 
I972I: iLavallev et'al. 1992) . We find that the two 
distributions are inconsistent at the ^92-98% con- 
fidence level. In order to assess the effect of uncer- 
tainties on the orbital periods of the M31 systems 
we repeated the test several times by randomly 
sampling periods assuming that they follow a Nor- 
mal distribution with the mean and standard de- 
viation calculated from the values in Table 2. In 
all cases the results give very similar confidence 
levels. Our sample of M31 transients differs from 
the MW transients at the ^2a level; increasing 
the M31 transient population should amplify this 
result. 

3.4-2. Trends in period and B magnitude limit vs 
distance from the M31 nucleus 

Voss fc Gilfanovl (|2007l ) discovered an enhanced 
population of X-ray binaries within ^100" of the 
M31 nucleus. The radial distribution of "primor- 
dial", i.e. expected, LMXBs follo ws the distribu - 
tion of stellar mass in the galaxy (|Gilfanovll2004 ): 
however, the surplus population appears to fol- 
low the p1 distri bution expected of dynamically 



formed binaries ( Fabian et al.l Il975l ). Dynami- 



cally formed LMXBs are normal ly very rare out- 
side g lobular clusters; however, IVoss &: Gilfanovl 
(|2007l ) propose that the bulge of M31 is sufficiently 



large and dense to allow the formation of a sig- 
nificant number of dynamically formed binaries. 
However, the stellar velocities in the M31 bulge 
are ^ 5-10 times higher than in globular clusters; 
hence IVoss fc Gilfanovl (|2007n expect only short 
period binaries to survive, with most of the sur- 
vivors containing black holes. 

Since the number of dynamically formed black 
binaries is strongly dependent on the distance 



from the center (rgc), one might expect to observe 
some correlation between orbital period and ra- 
dial distance for our transients. In the top panel 
of Fig. [5] we plot rgc vs. P for each transient, ex- 
cluding T8. The bottom panel of Fig. [S] shows 
F435W (~B band) detection limit for each tran- 
sient; T7 lies in a diagonally-shaded region be- 
cause it showed no X-ray variability between the 
HST on and off observations, hence the B mag- 
nitude was limited by the brightest star within 
the 3(7 location ellipse rather than the local back- 
ground. The sensitivity increased steadily with 
increasing rgc, varying by 2.5 magnitudes (i.e. a 
factor 10). 

We employed the Cox hazard test and the 
Kendal tau test, both implemented in ASURV to 
deal with bivariate data, to compare the popula- 
tion located 20" < rgc <150" with the popula- 
tion at 150" < Tgc <350"; the probability that 
these two populations were drawn from the same 
population (po) is ^0.12. Limiting rgc to 300", 
thereby removing T2, and re-running the Cox haz- 
ard and Kendal tau tests produces the same result. 
By contrast, extending r-gc to 400" for the sec- 
ond group, thereby including Tl, eradicates any 
trace of a correlation. Therefore, a more substan- 
tial transient population is required to confirm this 
suspected correlation. 

3.4.3. Decay profile vs. X-ray luminosity and es- 
timated period 

X-ray transients are expected to exhibit ex- 
ponential or linear decay, depending on whether 
or not the disc is fully ionized; this depends 
on the orbital period and peak X-ray luminosity 
(|Shahbaz et al.lll998l) . Hence we used Fig. 1 of 



that paper to determine whether the decays of our 
targets should be exponential or linear. 

With two exceptions, all of our targets are ex- 
pected to exhibit exponential decay; T7 is con- 
sistent with either linear or exponential decay, 
while Til s hould decay linearly according to 
(|Shahbaz et al. 1998 ). 

Interestingly the decay of Til was definitely 
more exponential (x^/dof = 2/5) than linear 
(x^/dof = 93/5); we suggest that the true peak lu- 
minosity of Til was considerably higher than the 
observed peak 0.3-10 keV peak of '-^4x10'^^ erg 
s~^. The true orbital period of Til could then be 



significantly shorter than our estimate of 150±30 
hr. 

Furthermore, the T5 outburst appears to have 
decayed linearly, even though exponential decay is 
expected; no observations were made in the ~170 
days prior to the observed outburst, hence we may 
have caught the linear tail of a longer and brighter 
outburst. A brighter peak X-ray luminosity for T5 
would mean a shorter period; however, this would 
not affect our results, since we only have an upper 
limit for T5. 

4. Summary and conclusions 

We report new orbital period estimates for 
12 X-ray transients observed by Chandra and 
HST. We have identified counterparts to 4 X- 
ray sources, 2 of those are new to this work; we 
have estimated 4cr upper limits to the B mag- 
nitudes of the remaining 8 systems. Period es- 
timates were previously published for 5 sources 
(Tl, T3, T5, TI P, T12) , using the VP M9 4 em- 
pirical relation (IWilliams et all l2Q05all3 E l2006t 
iBarnard et aP l201l[ ). and Noorae et al. (2012) 
have already applied the Ml 2 relation to T9. We 
obtained periods that ranged from <0.17 hours to 
<175 hours for black hole LMXBs using the M12 
relation, and <0.4-490±90 hours for neutron star 
LMXBs using the VPM94 relation. 

The M31 transients that we have followed up 
appear to favor shorter periods than the Galactic 
black hole transients; if all of our transients con- 
tain BHs, then the probability for the M31 and 
MW binaries being drawn from the same popula- 
tion is <0.08. This result may be explained by an 
excess of X-ray binaries dynamically formed in the 
high density bulge of M3 1; such systems requir e 
short periods to survive (jVoss &: Gilfanovl 120071) . 
However, some of the transients with the shortest 
periods may contain neutron stars. 

We searched for a correlation between the pe- 
riod and distance from the nucleus that would 



furthe r support the hypothesis of lVoss fc Gilfanov 
(J2007l ). We found a hint of a correlation, at the 
^90% confidence level. 

Further observations of M31 transients with 
Chandra and HST will add significance to these 
relations. 
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Fig. 1. — Locations of the 12 transients, superposed on a B band image from the LGS (Massey et al., 2006). 
The area enclosed by the circle is within 100 " of the nucleus: this reg ion is thought to produce X-ray binaries 
at an enhanced rate, due to its high density IVoss fc Gilfanovl (|2007h . 
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Fig. 3. — Difference images for the 12 transients, obtained by subtracting the HST off image from the HST 

on image. Each image is 5"x5". White is positive, and black is negative. The dotted ehipses indicate 

la uncertainties, while the solid ellipses indicate Scr uncertainties in the location of the X-ray source. We 

subtracted a V-Band off image from a B-band on image for Tl; there are many black stars because they are 

brighter in V than in B. 
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Table 1: Journal of HST ACS observations. For each observation we provide the HST ID, date, exposure 
time in seconds, RA and Dec of pointing, and target. We also show the closest Chandra ACIS observation, 
and the time of the Chandra observation relative to the HST observation. 



Obs 


ID 


Date 


Exp 


RA 


Dec 


Target 


Chan Obs 


Sep. 


1 


8vp02010 


2003-12-25 


2200 


00:42:56.03 


-^41:12:19.0 


TlOon 


4680 


+2 days 


2 


8vp03010 


2004-01-23 


2200 


00:42:43.86 


+41:16:30.1 


T6 on 


4681 


+8 days 


3 


8vp04010 


2004-06-14 


2200 


00:42:33.39 


+41:17:43.4 


T3 on 


4682 


-22 days 


4 


8vp05010 


2004-08-15 


2200 


00:42:41.50 


+41:17:00.0 


T5 on 


4719 


+29 days 


5 


8vp07010 


2004-10-02 


2200 


00:42:56.04 


+41:12:19.1 


TlOoff 


4721 


+2 days 


6 


8vp08010 


2004-11-02 


2200 


00:43:09.93 


+41:23:32.6 


T12 on 


4721 


-29 days 


7 


8vp09010 


2004-11-22 


2200 


00:42:41.49 


+41:17:05.2 


T5 off 


4723 


+13 days 


8 


8vpl0010 


2005-01-01 


2200 


00:43:09.91 


+41:23:32.5 


T12off 


4723 


-27 days 


9 


j9ju01010 


2006-02-10 


4360 


00:42:52.38 


+41:16:49.2 


T8 on 


7136 


—35 days 


10 


j9ju02010 


2006-07-09 


4360 


00:42:33.89 


+41:23:31.1 


T4 on 


7139 


+22 days 


11 


j9ju03010 


2006-07-15 


4352 


00:42:51.17 


+41:19:17.9 


T7 on 


7139 


— 16 days 


12 


j9ju04010 


2006-08-27 


4344 


00:42:16.08 


+41:19:26.1 


Tl on 


7139 


-27 days 


13 


j9ju05010 


2006-10-25 


2148 


00:42:17.34 


+41:15:37.0 


T2on 


7140 


-32 days 


14 


j9ju06010 


2007-01-10 


4672 


00:42:52.38 


+41:16:49.2 


T8off 


8183 


+4 days 


15 


j9ju07010 


2007-01-12 


4672 


00:42:33.88 


+41:23:31.2 


T4off 


8183 


+2 days 


16 


j9ju08010 


2007-01-12 


4672 


00:42:51.14 


+41:19:18.0 


T7ofr 


8183 


+2 days 


17 


9udl7010 


2009-08-25 


4360 


00:42:17.31 


+41:15:37.5 


T2off 


10555 


—55 days 


18 


b9dl5010 


2010-01-21 


4360 


00:42:53.15 


+41:14:22.9 


T9 on 


11278 


+14 days 


19 


b9d20010 


2010-07-20 


4360 


00:42:53.15 


+41:14:22.9 


T9ofr 


11840 


—3 days 


20 


jbdlOlOlO 


2010-12-09 


5472 


00:42:57.00 


+41:20:05.0 


Til on 


12162 


-1 day 


21 


j9ju08010 


2011-07-20 


4672 


00:42:51.17 


+41:10:17.9 


Til off 


12972 


-5 days 
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Table 2: For each object, we provide the X-ray position, with la uncertainties, along with the 0.3-10 kcV 
luminosity during the closest Chandra observation to the HST on observation, and the estimated 2-10 kcV 
luminosity at the time of the HST on observation; the luminosities are normalized to 10^^ erg s~^. We 
then provide the B band apparent Vega magnitude, and the derived absolute V band magnitude. We finally 
estimate the orbital period using the VPM94 method (Pns)j a-nd M12 method (Pbh)- 



Targ 


RA 


ctra 


Dec 


CTDoc 


-^Ob 


-^Est 


B 


My 


Pm 


Pbh 


1 


00:42:16.106 


0.015 


+41:19:26.36 


0.27 


8.9±2.3 


0.3 


>28.7 


2.45 


<10.6 


<11.6 


2 


00:42:17.330 


0.012 


+41:15:37.48 


0.11 


12.8+1.2 


2.89 


24.7+0.3 


-0.2 


120+5 


mtl 


3 


00:42:33.426 


0.012 


+41:17:03.56 


0.14 


3.7+0.3 


1.59 


>27.2 


2.3 


<3.6 


<4.5 


4 


00:42:33.883 


0.016 


+41:23:31.49 


0.19 


2+0.2 


0.88 


>28.9 


4 


<0.4 


<0.97 


5 


00:42:41.805 


0.010 


+41:16:35.94 


0.10 


5.8+0.8 


2.52 


>26.6 


1.7 


<6.0 


<5.6 


6 


00:42:43.834 


0.011 


+41:16:30.60 


0.11 


0.94+0.12 


0.53 


>26.3 


1.4 


<35 


<31 


7 


00:42:51.199 


0.012 


+41:19:17.81 


0.13 


0.51+0.1 


0.35 


>25 


2.9 


<380 


<175 


8 


00:42:52.433 


0.010 


+41:16:48.75 


0.10 


0.35+0.08 


0.29 


>23.9 


-1 


<2400 


<1000 


9 


00:42:53.181 


0.010 


+41:14:22.83 


0.10 


38+9 


18.62 


23.91+0.08 


-0.99 


61+11 


1811 


10 


00:42:56.049 


0.012 


+41:12:18.64 


0.13 


3+0.2 


1.77 


24.52+0.02 


-0.38 


163+8 


911 


11 


00:42:56.991 


0.017 


+41:20:05.35 


0.31 


0.39+0.10 


0.2 


24.87+0.09 


-0.03 


490+90 


uolf. 


12 


00:43:09.963 


0.030 


+41:23:32.56 


0.43 


1.65+0.12 


0.53 


> 27.9 


3 


<2.8 


<4.8 
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Table 3: Spectral fit parameters for transients with >200 net source counts in their spectra, assuming an 
absorbed disk blackbody model. For each source we give the net source counts, absorption, inner disk 
temperature and x^/dof for the best fits. We also provide the ratio of luminosities in the 2-10 and 0.3-10 
keV bands. We note that the T9 spectrum is piled up; this is accounted for in Noorae et al. (2012) 

Object Net source counts A^h / 10^^ atom cm"^ kT / kcV x^/dof i2-io/i0.3- 10 



Tl 


374 


T2 


756 


T5 


330 


T9 


1135 



2.5±1.9 


1.5±0.3 


10/15 


0.55 


0.67 


0.54±0.04 


35/31 


0.25 


0.67 


1.1±0.2 


4/12 


0.59 


0.67 


0.8 


64/49 


37 
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